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electromagnetic field generated by a radio base station 
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Description 

[0001] This invention refers to the means for planning 
telecommunication systems on radio carriers, and in 
particular refers to a method for determining the values 
of an electromagnetic field generated by a radio base 
station in an urban environment. 
[0002] The precise evaluation of electromagnetic field 
levels near radio base stations in order to check that 
legal restrictions have been observed, is today one of 
the fundamental aspects in the designing of mobile radio 
networks. Since the on-site measurement of electro- 
magnetic field levels is often prohibitive cost-wise, and 
furthermore is not very selective when determining irra- 
diating contributions, software instruments are becom- 
ing ever more frequently used to estimate the field lev- 
els. 

[0003] Unfortunately the complexity of urban environ- 
ments makes it difficult to solve the problem with theo- 
retical calculations based on Maxwell's equations. 
[0004] For this reason, and due also to the approxi- 
mate description of the buildings of the environment un- 
der examination, methods have been used that estimate 
the electromagnetic field as the sum of contributions of 
reflected and diffracted rays that propagate in a straight 
line from the source. 

[0005] An electromagnetic field evaluation applica- 
tion, as for example that described in "The Mobile Radio 
Propagation Channel" 2nd edition, by J. D. Parsone, Ed. 
John Wiley & Sons LTD, generally works in the following 
way: starting from a data base containing a vector map 
of the building and a data base of the irradiating sources, 
it calculates the optical paths of the signal that reaches 
a generic area of the urban environment. Subsequently, 
the optical paths calculated with the numerical calcula- 
tion models can be used to find the values of the elec- 
tromagnetic field, which can then be used to estimate 
the environmental impact that is a fundamental element 
in the design of a cellular network. 
[0006] In the urban environment, the calculation of the 
optical paths is the most time-consuming part of the ap- 
plication. The complexity of the algorithm, in fact, is such 
that even the analysis of small areas of a few hundred 
metres takes far too long using present-day calculation 
instruments, unless optimisation procedures are used. 
For this reason, various techniques have been devel- 
oped to introduce approximations that make it possible 
to perform this type of calculation in a reasonable 
amount of time. These techniques can be subdivided 
into two different categories: 

forward ray-tracer algorithms; 
backward ray-tracer algorithms. 

[0007] In forward ray-tracer algorithms, a finite 
number of rays that are isotropically irradiated from the 
transmitter and their propagation into space is consid- 
ered, taking into consideration possible reflections and 



diffractions with elements in the environment. To estab- 
lish whether a ray reaches the receiver, the receiver is 
taken to be located at the centre of a sphere constituting 
the capture volume, and it is determined whether the ray 

5 intersects the sphere. This type of algorithm is suitable 
for studying propagation in a two-dimensional environ- 
ment, whereas it is not very efficient when studying a 
three-dimensional environment. This is due to the con- 
siderable number of rays that must be irradiated isotrop- 

io ically into the solid angle, regardless of the direction in 
which there are objects that can interact. An operation 
that involves a great number of unnecessary calcula- 
tions. 

[0008] On the other hand, the backward ray-tracer al- 
15 gorithms calculate which rays connect a transmission 
source with a reception point, taking into account that 
they may, during their path, be reflected and diffracted 
a considerable number of times. 
[0009] The most critical aspect of a ray-tracer algo- 
20 rithm, from the point of view of calculation time, is the 
"visibility test" performed between source and receiver, 
i.e. the check that there are no obstacles between the 
source and the receiver that may interrupt the rays. To 
perform this test, the most efficient algorithms are ar- 
25 ranged to carry out two main operations: 

generate a structure, called "visibility tree", in which 
the visibility relationships between the various ob- 
jects under examination, are stored; 
30 - calculate optical paths. 

[0010] The first operation builds up the visibility rela- 
tionships between the various objects in the urban en- 
vironment under examination. In practice, it defines the 
objects that can be reached by a ray that is reflected or 
diffracted by another object. The building of the visibility 
tree involves the problem of generating two-dimensional 
(2-D) images of three-dimensional (3-D) environments, 
which is called "image synthesis", and is performed by 
electronic processors. In actual fact, these techniques 
are very refined and only the parts required for the vis- 
ibility test need be extrapolated forthe propagation mod- 
els, as will be explained later. 

[0011] The second operation of a ray-tracer algorithm 
consists in searching for all the possible optical paths 
between transmitter and receiver, using the visibility 
tree. 

[0012] The building of the visibility tree, which is the 
object of the first operation, has an algorithmic complex- 
ity that grows exponentially with the number of reflec- 
tions requested. For this reason, as already mentioned, 
electronic image processing techniques are used to 
speed up the building and reduce the complexity of the 
tree. 

[0013] The number of possible surfaces interacting 
with the rays irradiated from the source can be reduced 
by evaluating the angle formed between the segment 
that joins the source to one of the surface points under 
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examination and the surface point's perpendicular, ac- 
cording to the "back-face culling" technique, described 
on pages 663-664 of the book entitled "Computer 
Graphics: Principles and Practice" 2 nd edition, by J. D. 
Foley et al. When this angle exceeds ninety degrees the 
surface can be ignored in that it is physically impossible 
for a ray irradiated from the source to be reflected on it. 
[0014] The number of surfaces can be further reduced 
by identifying those surfaces that are completely hidden 
by other objects. To do this, the article "Efficient ray-trac- 
ing technique for three-dimensional analyses of propa- 
gation in mobile communications: application to picocell 
and microcell scenarios," by M.F. Cattedra et al., IEEE 
Antennas & Propagation Magazine, vol. 40 pp. 15-28 
April 1 998, presents a series of geometric relations to 
verify whether a surface completely blacks out another 
surface with respect to the source. The disadvantage of 
this technique resides in the fact that it is necessary to 
define a considerable number of angular relationships 
that are a function of how the two surfaces are posi- 
tioned in space (vertical-vertical, vertical-horizontal, 
vertical-slanting, etc.), and in the fact that all special cas- 
es (surfaces seen under an angle which comprises the 
direction 0 or 2n) must be dealt with. 
[0015] To overcome this problem, F. Brunello, D. Dis- 
co and D. Gambin, in "An acceleration technique using 
a 3D representation for ray tracer in a urban environ- 
ment', IEEE Antenna and Propagation Symposium 
2000, suggest that a perspective representation centred 
on the source be used and implemented with well- 
known computer graphics processing techniques. 
[0016] By representing the surfaces starting with 
those furthest away from the source, the closest surfac- 
es are superimposed on the others making them disap- 
pear from the overall image. The "surviving" surfaces 
are the ones that will be considered by the ray tracer. 
[0017] The main problem with this technique is the in- 
troduction of segment distortion into the perspective 
representation, if the segments are represented through 
the segment joining the ends: non-linear transforma- 
tions are in fact used that consequently distort the sur- 
faces. 

[0018] The considerable use of trigonometric func- 
tions also makes the relative processing particularly dif- 
ficult from a computation point of view. 
[0019] The method for determining the values of an 
electromagnetic field generated by a radio base station 
in an urban environment, which is the subject of this in- 
vention, eliminates the aforesaid disadvantages and 
solves the technical problems described. It gives an ac- 
curate prediction, while the processing time required by 
the data processing equipment to make the calculation 
is much less. This is due a) to the reduction in the struc- 
ture that stores the visibility relationships between the 
various objects under examination, i.e. the visibility tree, 
and b) to the simplification of the visibility test used to 
determine the optical paths. 

[0020] The subject of this invention is a method for 



determining the values of an electromagnetic field gen- 
erated by a radio base station in an urban environment, 
as described in the characterising part of claim 1 . 
[0021] The foregoing and other characteristics of this 
5 invention will be made clearer by the following descrip- 
tion of a preferred form of the invention, given byway of 
non-limiting example, and by the annexed drawings in 
which: 

10 - Fig. 1 is a schematic representation of an urban en- 
vironment with a transmitter, a receiver and several 
buildings; 

Fig. 2 is a schematic representation of an urban en- 
vironment that illustrates the construction of the 
15 transmitter images; 

Fig. 3 gives an example of a visibility tree; 
Fig. 4 illustrates a projection according to a tech- 
nique called "z-buffer"; 

Fig. 5 shows an example of matrixes F and 2 ac- 
20 cording to the z-buffer technique; 

Fig. 6 is the flow diagram that illustrates the method 
of this invention. 

[0022] The evaluation method proposed herein uses 

25 backward ray-tracer algorithms, which considerthe rays 
that connect a transmission source T with a reception 
point R, as illustrated in Fig. 1. In this case, the direct 
path from T to R is not possible as building C is in the 
way, while the paths exploiting the reflection on the sur- 

30 faces of buildings B, C and D and the diffraction on the 
edge of building A are possible. 
[0023] The environment illustrated in Fig. 2 contains 
the transmission source Tx and the reception point Rx 
that are surrounded by four surfaces W.,, W 2 , W 3 and 

35 w 4 . Starting from Tx, image I^W^ is constructed with 
respect to the surface then the image l 2 (W 3 ) of im- 
age li(W.,) with respect to surface W 3 , then image i 3 
(W 4 ) of image l 2 (W 3 ) with respect to surface W 4 and fi- 
nally image l 2 (W 2 ) of image l 2 (W 3 ) with respect to sur- 

40 face W 2 . Starting from reception point Rx, the rays are 
traced to the last image until they reach the point of in- 
tersection with surface W 2 , thus obtaining the reflection 
point P 1t then from P 1 to the penultimate image l-jfW^ 
to obtain P 2 , and so on until source Tx is reached. 

45 [0024] Whenever a ray is traced between Rx and a 
reflection point or between two reflection points, it is nec- 
essary to check whether there are any obstacles in the 
way that may interrupt the ray, i.e. the visibility test must 
be carried out, which as already mentioned involves a 

so considerable number of calculations if all the possible 
obstacles are to be taken into consideration. 
[0025] To run this test, a visibility tree must be built. 
Fig. 3 shows a typical visibility tree structure between 
Tx and Rx. To make comprehension easier, the tree has 

55 been developed as far as the third reflection only for the 
paths illustrated in the Fig, 2 environment. To be more 
precise, images I^W,,) are associated with the first re- 
flection, images l 2 (W n ) are associated with the second 
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reflection and images l 3 (W n ) are associated with the 
third reflection. Each node contains the image's coordi- 
nates and a reference of the surface that generated it. 
[0026] As has been stated previously, the visibility 
tree is rather complex; in fact the number of images to 
be calculated for n reflections is N • (N - 1) 5 , N being the 
number of objects in the environment under examina- 
tion. To speed up the building of this tree, the method 
referred to in the invention uses both the "back-face cull- 
ing" technique mentioned previously, and a technique 
for reducing the depth of the visibility tree that exploits 
the z-buffer projection method. 
[0027] The "back-face culling" technique is used to 
exclude all the objects that do not intersect the ray re- 
flected from the surface under consideration. To check 
that an object does not constitute an obstacle to the ray, 
simply check that at least one of the following conditions 
exists: 

the object is not in the semi-space delimited by the 
plane of the reflecting surface, outside the building; 
the reflecting surface is not in the semi-space de- 
limited by the plane of the object 1 s surface, outside 
the object itself. 

[0028] The technique can be applied even if the ray 
comes from a point source, and not from a reflecting sur- 
face. In this case, simply check that the source is not in 
the semi-space delimited by the plane of the objects 
surface. 

[0029] This technique does not determine which ob- 
jects in the environment constitute an obstacle to the 
rays, but those that do not. 

[0030] The projection technique, called "z-buffer", is 
based on the idea of perspectively projecting the objects 
in the environment on the plane of the reflecting surface 
W n that has generated image ! m (W n ) of the source, the 
visibility of which is to be studied, as illustrated in Fig. 
4. In this figure, object OD is projected on plane PP, on 
which the reflecting surface FQ lies, parallel to plane x, 
y of the reference system x, y, z, the origin of which is 
in image I. 

[0031] In terms of implementation, the z-buffer tech- 
nique includes two matrixes, Fand Z, with dimensions 
m • n, as shown in Fig. 5. 

[0032] Each cell in matrix F represents a geometric 
element of plane PP (Fig. 4) and when this element co- 
incides with a projection element of an object in the en- 
vironment, a univocal reference is stored, e.g. a pro- 
gressive whole number corresponding to the object that 
has generated the projection. 

[0033] Similarly, each cell in matrix Zcontains coordi- 
nate z of the object's corresponding element, the pro- 
jection of which is represented in matrix F. 
[0034] To be more precise, if a polygonal object is to 
be represented on matrix F, its vertices aae projected 
onto plane PP thus obtaining the respective cartesian 
coordinates u, v for each one. The number indicating 



the polygon is entered in cells f uv of F; this number is 
also entered in all the cells that join the f uv cells relating 
to the vertices, placed on straight segments, and in eve- 
ry cell contained in the perimeter that has been just 
5 traced so as to complete the representation of the pol- 
ygon on the matrix. 

[0035] The projections of several objects in the envi- 
ronment may fall on the same geometric elements of the 
projection plane PP, but only the z coordinates of the 

10 object nearest the projection plane must be found in the 
corresponding cells in the Zmatrix. Consequently, after 
an object has been perspectively projected on the plane 
that has generated the image, the visibility of which is 
being studied, the z coordinate of the point of the object 

15 generating it is calculated for each geometric element 
in the projection. If the value is less than that already 
possibly in the cell as a result of a previous projection, 
then the new z value is entered, and the number indi- 
cating the object is entered in the corresponding cell of 

20 the F matrix. 

[0036] Finally, in order to take into account all the ob- 
jects seen in the image under examination, and only 
those objects, all the cells relating to geometric ele- 
ments that do not belong to the FQ surface area that 

25 generated the projection plane PP must be eliminated 
from the matrixes. 

[0037] In the method referred to in the invention, z- 
buffer projection is used to build the visibility tree. To do 
this it is necessary to make the hypothesis that all the 
30 objects in the urban environment can be estimated by 
joining flat surfaces. 

[0038] If this hypothesis applies, then the buildings 
can be described with a set of elementary polygons (tri- 
angles, rectangles). 
35 [0039] In this way, all the objects that can be reached 
by a ray reflected from an FQ polygon coincide with 
those in the projection plane for a system of coordinates 
in which: 

40 a) The origin of the coordinated axes system coin- 
cides with image I of the source with respect to the 
plane on which the FQ polygon lies, 
b) The projection plane Z is equal to the distance of 
the FQ polygon from image I. 

45 

[0040] The construction of the visibility tree is there- 
fore reduced to three basic steps: 

Back-face culling to determine which surfaces face 
so the reflection plane. 

Rotation and translation of the axes system to sat- 
isfy conditions a) and b). 
z-buffer projection. 

55 [0041] The z-buffer technique is used in a new way to 
reduce the visibility tree depth by one level. By evaluat- 
ing the distances in matrix Z it is possible to identify the 
objects nearest the reflection plane. Consequently a re- 
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ception point R is "visible" if its z coordinate is less than 
the z coordinate relating to the geometric element on 
which the projection of R falls. 

[0042] The evaluation method that uses the above 
techniques is shown in the flow diagram Fig. 6. 
[0043] The method starts at step 1 , then goes to step 
2, where the urban environment under examination is 
loaded by creating a list of object surfaces, each one 
identified by a progressive number, by the coordinates 
of its vertices, and by a matrix for the rotation and trans- 
lation operations in the visibility calculation, etc. 
[0044] The next step, number 3, is where the building 
of the visibility tree starts, the first operation of which, 
step 4, is the back-face culling of each element in the 
surfaces list, to initially exclude any objects that are def- 
initely not visible. 

[0045] Step 5 creates the list of source images and 
enters this list in the visibility tree at level 1 . 
[0046] Step 6 initialises the number of reflections, set- 
ting it to 1 . 

[0047] Step 7 checks to see whether the number of 
reflections is less than a maximum pre-established 
number of reflections, and depending on the result, the 
procedure may go in different directions. 
[0048] In particular, if it is true that the maximum 
number of reflections has not been reached, answer 
"yes" to the check, then variable J ranging from 1 to the 
number of images for each level of the visibility tree, is 
initialised at 1 , step 8. 

[0049] Step 9 checks to see whether variable J is less 
than or equal to the number of images of the current 
level of the visibility tree, and depending on the result 
the procedure may go in different directions. 
[0050] In particular, if it is not true that J is less than 
or equal to the number of images, answer "no" to the 
check, then the number of reflections is increased by 
one unit, step 1 0, and then the procedure continues to 
the next visibility tree level, step 1 1 , to then return to step 
7. 

[0051] If, however, it is true that J is less than or equal 
to the number of images, answer "yes" to the check, 
then the operations included in the back-face culling 
technique are performed for each surface of the objects 
in the list, step 12, thus obtaining a short list containing 
only the surfaces that can be seen from the image. 
[0052] Step 1 3 initialises variable K at 1 , which is used 
atthe next step, 14, to examine all the surfaces just iden- 
tified, and checks to see whetherthe number is less than 
or equal to the number of surfaces, and depending on 
the result the procedure may go in different directions. 
[0053] In particular, if it is not true that K is less than 
or equal to the number of surfaces, answer "no" to the 
check, then variable J, representing the number of im- 
ages, is increased by one unit, step 15, to then return to 
step 9. 

[0054] If, however, it is true that K is less than or equal 
to the number of surfaces, answer "yes" to the check, 
step 1 6 is carried out to translate the surface to make it 



coincide with the origin of the reference system, and 
step 1 7 is performed to rotate the same surface to make 
it parallel with the x, y plane. 

[0055] For each translated and rotated surface, the z- 
5 buffer technique performs projection, step 18, and then 
the visibility test of the reception points is run using their 
projection on the z-buffer plane and comparison with the 
z coordinate, step 19. 

[0056] Atthe next step, step 20, it is therefore possible 

10 to create the first section of the path, starting from the 
reception point that is shown to be visible in the previous 
visibility test, up to the first level reflection point. 
[0057] Variable K is then increased by one unit in or- 
der to examine the next surface, step 21, and then it 

'5 returns to step 14. 

[0058] When step 7 shows that the number of reflec- 
tions examined has reached the maximum pre-estab- 
lished number, answer "no" to the check, the procedure 
moves onto the direct visibility test between transmitter 

20 and receivers, step 22, and then the visibility test to com- 
plete the optical paths, the first section of which has al- 
ready been calculated, step 23. 
[0059] The method ends at step 24. 
[0060] As can be seen, the z-buffer graphics acceler- 

25 ation technique has not only be used to build the visibility 
tree (step 18), but also to calculate the first part of the 
optical paths (step 20), eliminating the visibility compar- 
ison between the reception point and the last point of 
reflection. By doing so the depth of the visibility tree is 

30 reduced. 

[0061] Obviously this description is given as a non- 
limiting example. Variants and modifications are possi- 
ble, without emerging from the protection field of the 
claims. 

35 

Claims 

1. Method for determining the values of an electro- 
40 magnetic field generated by a radio base station in 
an urban environment, using a "backward ray-trac- 
er" algorithm to calculate which rays connect a 
transmission source (T) with a reception point (R), 
taking into consideration that the rays can be re- 
45 fleeted and diffracted by the surface of objects in 
the environment during their path, characterised 
by the fact that the "back-face culling" technique is 
used to exclude all the objects that fail to intersect 
the rays reflected from the surfaces considered, 
so and the "z-buffer" projection technique is used to 
reduce by one level the structure, called the "visi- 
bility tree", of the visibility relationships between the 
various objects under examination. 

55 2. Method as in claim 1 , characterised by the fact that 
it includes the following steps: 

start (step 1); 
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loading of the urban environment by creating a 
list of object surfaces, each one identified with 
a progressive number, its vertex coordinates, 
and a matrix for rotation and translation opera- 
tions (step 2); 5 
start of visibility tree construction (step 3); 
back-face culling for each element in the list of 
surfaces, in order to exclude the objects that 
are definitely not visible (step 4); 
creating the list of source (T) images and its en- 10 
try at the first level of the visibility tree (step 5); 
initialising the number of reflections (step 6); 
checking to see whether the number of reflec- 
tions is less than the maximum pre-established 
number of reflections (step 7); is 
if the maximum number of reflections has not 
yet been reached, the first variable (J), which 
varies from 1 to the number of images for each 
visibility tree level, is initialised (step 8); 
checking to see if the first variable (J) is less 20 
than or equal to the number of images of the 
current visibility tree level (step 9); 
if the first variable (J) is not less than or equal 
to the number of images, the number of reflec- 
tions is increased by one unit (step 1 0); 25 
moving on to the next visibility tree level (step 
11), and then checking to see whether the 
number of reflections is less than the maximum 
pre-established number of reflections (step 7); 
if, however, the first variable (J) is less than or so 
equal to the number of images, the operations 
envisaged in the back-face culling technique 
are carried out for each surface of the objects 
in the list (step 1 2), thus obtaining a list contain- 
ing only the surfaces that can be seen from the 35 
image; 

initialising of the second variable (K), (step 13); 
checking to see whether the second variable 
(K) is less than or equal to the number of sur- 
faces identified (step 14); 40 
if the second variable (K) is not less than or 
equal to the number of surfaces, the first vari- 
able (J) (step 15) is increased by one unit, and 
then the first variable (J) is checked again to 
see whether it is less than or equal to the 45 
number of the images of the current visibility 
tree level (step 9); 

if the second variable (K) is less than or equal 
to the number of surfaces, the surface is trans- 
lated in order to make the image coincide with so 
the origin of a reference system (x, y, z) (step 
16); 

rotating of the same surface in order to make it 
parallel to plane (x, y) of said reference system 
(step 17); " ss 

for each translated and rotated surface, projec- 
tion is made according to the z-buffer tech- 
nique, (step 18); 



visibility test of the reception points (R) by pro- 
jecting them onto a z-buffer plane (PP) and 
comparing the (z) coordinate (step 19); 
creating the first section of the ray's path, start- 
ing from the reception point (R) that is "visible" 
in the previous test, to the first level reflection 
point (step 20); 

increasing of the second variable (K) by one 
unit so as to examine the next surface (step 21), 
and then checking again to see if the second 
variable is less than or equal to the number of 
surfaces identified (step 14); 
if the number of reflections examined has 
reached the maximum pre-established number 
(step 7), the direct visibility test between the 
source (T) and the receiver points (R) is carried 
out (step 22); 

running of the visibility test to complete the op- 
tical paths, the first section of which has already 
been calculated (step 23); 
end (step 24). 

3. Method as in claim 1 or2, characterised by the fact 
that, in order to carry out the visibility test, the visi- 
bility tree is built between source (T) and reception 
point (R), in which the source is connected to an 
initial set of nodes representing the images (I^W,,)) 
associated with the first reflection, each node of the 
first order being connected to several nodes of the 
second order representing the images (l 2 (W n )) of 
the second reflection and so on, each node contain- 
ing the coordinates of the image and a reference of 
the surface that generated it. 

4. Method as in any of claims 1 to 3, characterised 

by the fact that in order to project according to the 
z-buffer technique, objects in the environment (OD) 
are projected perspectively onto plane (PP) of the 
reflecting surface (FQ) that has generated the im- 
age (l m (W n )) of the source (T), the visibility of which 
is to be tested, said plane (PP) being parallel to the 
plane (x, y) of the (x, y, z) reference system, the or- 
igin of which is in the image (I). 

5. Method as in any of claims 1 to 4, characterised 
by the fact that a first matrix (f) and a second matrix 
(Z) of dimensions m • n are built according to the z- 
buffertechnique, each cell of thefirst matrix (F), rep- 
resenting a geometric element of the plane (PP) of 
the reflecting surface (FQ), containing a univocal 
reference to the object that generated the projec- 
tion, when this element coincides with an element 
of the projection of an object in the environment. 

6. Method as in claim 5, characterised by the fact that 
each cell of the second matrix (2) contains the co- 
ordinate (z) of the corresponding element of the ob- 
ject, the projection of which is represented in the 
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first matrix (F), 

7. Method as in claim 6, characterised by the fact 
that, if the projections (OP) of more than one envi- 
ronment object (OD) fail on the same geometric el- 5 
ements of the projection plane (PP), only the minor 

(/) coordinates, which correspond to those of the 
object nearest the projection plane (PP), are en- 
tered in the second matrix (2) cells, and the univocal 
object reference is entered in the corresponding cell io 
in the first matrix (F). 

8. Method as in claim 4 or any of claims 5 to 8 if re- 
ferred to claim 4, characterised by the fact that all 

the cells relating to the geometric elements that do is 
not belong to the surface area (FQ) that generated 
the projection plane (PP) are eliminated from the 
matrixes (F, 2). 
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